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I. INTRODUCTION 

Electron scattering and transport through quantum 
dots (QDs) in a semiconductor nanostructureii^ii^ii^i^ 
have been intensively studied recently. The spin- 
dependent transport properties are of particular interest 
for its possible applications, e.g., the QD spin valves^, 
the quantum logic gates using coupled QDs, as well as 
the spin-dependent transport in single-electron devices^, 
etc.. In such systems, the electron-electron exchange po- 
tential and the electron spin states have been utilized 
and manipulated^i^iiSiiii. A thorough quantitative un- 
derstanding of spin-dependent properties due to electron- 
electron interaction is therefore important for a successful 
construction of these devices. This subject has been in- 
vestigated in different issues, such as the utilization of the 
electron-electron scattering in determining the electron 
entanglement dynamicai^, the study of spin-flip scatter- 
ing in double QDsl2. and the scattering through a region 
of nonuniform spin-orbit coupling which can form a spin 
polarized bearoi^. Theoretically the transport through 
QDs has been studied by different approaches such as 
transfer matrix, nonequilibrium Green's functions, ran- 
dom matrix theory, as well as those methods built on the 
Lippmann-Schwinger (L-S) equation. 

In this work, we develop a theoretical method to study 
electron scattering through a quantum dot (QD) of N- 
electrons embedded in a two-dimensional (2D) semicon- 
ductor system. We construct the scattering equations in- 
cluding electron-electron interaction to represent the pro- 
cess of a 2D free electron scattered by the QD. The gener- 
alized multichannel Lippmann-Schwinger equation o^^i^^ 
are solved for this system by using the method of contin- 
ued fractions (MCF). The MCF is an iterative method to 
solve the integro-differential L-S equations, initially de- 
veloped for three-dimensional electron-atom (molecule) 



scattering in atomic physicsii. We show that this 
method, combining with the partial wave expansion tech- 
nique, is of a rapid convergency for the present problem 
in a 2D semiconductor system and therefore is efficient to 
obtain the scattering cross-sections. As an example, we 
apply this method to a one-electron QD case and obtain 
scattering cross-sections resulting from both the singlet- 
and triplet-coupled continuum states of two electrons (in- 
cident and QD electrons) during the collision. The results 
show that the scattering processes can be very different 
for singlet and triplet spin states which are mainly origi- 
nated from the different exchange interactions. From the 
difference of the scattering amplitudes resulting from the 
singlet and triplet couplings, we determine the spin-flip 
scattering cross-sections which exhibit a maximum as a 
function of scattering angle and the incident electron en- 
ergy. In a multichannel scattering, we study the inelastic 
scattering process in which the incident electron is scat- 
tered by a lower energy state of the QD and leaves behind 
the QD in an excited state. As expected, such an inelas- 
tic scattering cross-section is found much smaller than 
the elastic one. 



This paper is organized as follows. In Sec. II we 
present the Hamiltonian of the system. In Sec. Ill, we 
describe our general theoretical approach and the one- 
electron QD case is given as an example. In Sec. IV 
we show our numerical results for the scattering through 
a one-electron QD within both the one-channel and the 
multichannel models. The conclusion is presented in Sec. 
V. Moreover, the method of continued fractions is briefly 
described in Appendix A. The 2D partial wave expan- 
sions used in the numerical solution of the L-S equations 
are presented in Appendix B. 
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II. HAMILTONIAN OF THE SYSTEM 

The system under investigation consists of an incident 
2D free electron and a quantum dot of N electrons em- 
bedded in a 2D system. The incident electron is scattered 
by both the QD potential and by the confined electrons 
inside the QD. The Schrodinger equation of the system 
is given by 



(1) 



where r represents collectively the spatial and spin co- 
ordinates of the N electrons localized in the QD and 
rw+i = {xn+ItUn+i) and a^+i denote the spatial and 
spin coordinates of the incident electron. The total en- 
ergy of the system is £i, where the subscript i represents 
a set of quantum numbers required to specify uniquely 
the initial quantum state of the system. Explicitly, the 
total Hamiltonian of the system can be written as 



H = Ho{rN+i) + Hqu{t) + Vint(i'i,r2, tat, tat+i) , 

(2) 

where Ho{rN+i) = -h \/%_^^/2m* + VQD(rAr+i), 
Hqi:>{t) is the Hamiltonian of the QD of N electrons, 
and Vint is the interaction potential between the incident 
electron at r^^i and the N electrons in the QD 



N 



Mnt(ri,r2, ...,rAr,r7v+i) = — X! 



1 



(3) 



where eg is the dielectric constant of the semiconduc- 
tor material and m* is the electron effective mass. The 
Hamiltonian for an unperturbed QD is given by 



N 



2m* 



^QD(ri) +-r 



T N 



1 



^5 - 



(4) 

where the first term in the rhs of Eq. Q describes N in- 
dependent electrons in the QD of confinement potential 
VQD(r) and the second term gives the Coulomb interac- 
tions among these electrons. The eigenenergy and eigen- 
function of this N-electron QD are denoted by £„ and 
respectively. They are determined by the following 
Schrodinger equation 



(5) 



with n = 0, 1, 2, 3... . The ground state of the N-electron 
QD is labeled by n = and the excited states by n > 1. 
The eigenstates of the QD can be obtained using, e.g.^ the 
restricted or unrestricted Hartree-Fock (HP) methods^S. 



III. SCATTERING EQUATIONS INCLUDING 
ELECTRON EXCHANGE INTERACTION 

In order to extract scattering properties of the sys- 
tem (QD + incident electron), we can write the total 



wave-function of the system as a superposition of the 
QD wave-function $" and the incident electron wave- 
function, 



(6) 



n=0 



where ipni describes the wave-functions of the incident 
(scattered) electron in the continuum states correspond- 
ing to a quantum transition from an initial state i to 
a final state n. The operator A warrants the antisym- 
metrization property between the QD electrons and the 
incident electron, defined by, 



, N+l 

* p—1 



(7) 



where Vn+i.p is the permutation operator which ex- 
changes the electrons at rjv+i and Yp. From Eqs. H]), 
([2|) and ([6]), we obtain 

oc 

£:,^|^($"V™)). (8) 



The total energy of the system (the incident electron -|- 
QD) £i is composed by two parts. The first part is the 
kinetic energy of the incident (scattering) electron and 
the second is the energy of the N-electron QD in a par- 

k^ 

ticular configuration, i.e., £i = — f -\-ei — — ^ for 
different eigenstates of the QD (i, n = 0, 1, 2, ...) or differ- 
ent scattering channels. These different channels appear 
because the incident electron can probably be scattered 
inelastically, leaving the QD in a different state from its 
initial. A projection of Eq. ([5]) onto a particular QD 
state 1$™) leads to the following scattering equation for 
the incident electron. 



(V^ + fc^) V'm«(r) = 14m(r)V'™(r) (9) 



2m* 



n=0 



for i, TO = 0, 1, 2, where r = r^+i and Kn„ = y,^'„ -I- 
with the static potential and the exchange 
potential due the nonlocal interaction, giving by 



N 



C'Jr) = yQD(r)<5,„„-f-^(<i> 
and 



-A|r-rj| 



|$"), (10) 



2 1,2 



C"„(r)^™(r) = (i/o(r) 



™)($"|^'(a>"v™)) 



2to* 



^ 1 
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respectively, where A' = J2p=i(^^)^^^ ^Vn+i.p- In 
Eq. ([TO|l we have introduced a screening e^'*'''""'' ' on the 
direct Coulomb potential for two reasons: (i) the the 
ionized impurities in the semiconductor nanostructure 
and/or the external electrodes screen the direct Coulomb 
potential and (ii) at the |r| ^ oo limit the scattering 
potential should decay faster than l/|r|. The screening 
length is given by A~^. Notice that we do not consider 
the screening on the exchange potential because this po- 
tential is non-zero inside the QD only. Inclusion of the 
screening on the exchange potential in Eq. (|lip is possi- 
ble but it will not affect much our results and complicates 
the numerical calculation. 

The scattering equation is a system of coupled integro- 
differential equations. The corresponding generalized L- 
S equation for such a multichannel scattering problem is 
given by 

1pmi{r) = (pi{r)Srm 

CO „ 

+ Y, /dr'G(")(k„,r,r')F„„(r')Vm(r'), (12) 

n=0 

for i, m = 0, 1, 2 . . . 

with an incident plane wave (pi{r) = e'*'' '' = e*''"'^ in the 
x-direction. The Green's function G^°-'(k, r, r') in the 
above equation is 

G(°)(k,r,r') = -^(»/4)//«(fc|r-r'|), (13) 

where Hq^'' is the usual zero order Hankel's function^i. 

At |r| — > cx) limit, the asymptotic form of Eq. (fT^ for 
the scattered wave-function in a 2D system is given by 

2m* / i g+ikmT 

IrHoo ^ V krn y/T 

(14) 

where fk^,ki{Q) is the scattering amplitude 

h^,km = -\\[l{\^ra\T{E)\k,) (15) 

with 

oc „ 

(k™|r(ii;)|k,) =^ / dr'e-'^''-'-V™„(r')VAu(r')- 

The momenta of the initial and final states of the incident 
(scattered) electron are ki and km, respectively, and 9 is 
the scattering angle between them. It is evident from 
Eq. p2| and its boundary condition Eq. (fT4l) that the 
different scattering channels are coupled to each other 
through the interaction potential Vmn- 

In the above procedure in dealing with the electron 
scattering through a QD, both the electron-electron ex- 
change and correlation interactions are presented. How- 
ever, it is difficult to include a complete correlation ef- 
fect in a practical calculation. For that, besides an exact 



solution for the N-electron QD, a full sum over all the 
intermediate states n in the scattering equation [Eq. ([9])] 
is needed, which is a formidable task in a self-consistent 
calculation. In an alternative way, the correlation effects 
can be considered by adding an effective correlation po- 
tential in the scattering equationi^. In the present work, 
we focus on the exchange effects in the scattering process 
and limit the sum over n to a few lowest energy levels of 
the QD. For this reason, we prefer to call the nonlocal 
interaction potential in Eq. (fTTj) as exchange po- 
tential though the correlation are partially included in a 
multichannel treatment. 

The differential cross-section (DCS) for a scattering 
from initial state i {i.e. the incident electron of kinetic 

energy Ei — and the QD in the state Si) to final 

state m (i.e. Em — 2rr? ^^'^ ™ state m) is 
given by 

= (16) 

The integral cross-section (ICS) which is an energy de- 
pendent quantity can be found by 

^m^{E^)= f a,m{0)de. (17) 

Jq 

When the state of the QD remains the same i.e. [m — i), 
before and after the scattering, the process is called 
elastic. Otherwise, the scattering is inelastic. A pos- 
sible scattering is the so-called super-elastic scattering 
{E,n > Ei) where the incident electron is scattered out 
with a higher energy by an QD initially in an excited 
state. Because the different scattering channels are cou- 
pled to each other, we have to solve the multichannel L-S 
equation to obtain the scattering probabilities through 
different channels simultaneously for the same total en- 
ergy of the system. 

There are different numerical methods to solve the 
above coupled integro-differential L-S equations. In this 
work, we use the so-called method of continued fractions 
(MCF, see Appendix A) which was originally developed 
in three-dimensional formulation for electron-ato m^^i 
and electron-molecule^ scatterings at the single- and 
multi-channel level of approximations. Here, we ap- 
ply this method to electron-QD scattering in a two- 
dimensional semiconductor system. The MCF is an it- 
erative method to solve the L-S equation. The advan- 
tage of this method lies on its rapid convergency and 
its unnecessity of a basis function for expansion of the 
continuum wave-functions. Using the MCF, we can ob- 
tain the T matrix and consequently the DCS according 
to Eqs. (|15p and p6)) . The two-dimensional integrations 
on the interaction potentials in Eqs. (jlOp . (Ilip . and (fT^ 
are simplified by using partial wave expansion which is 
described in Appendix B. 
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A. One-channel approximation 



with 



When a quantum dot is initiaUy in its ground state 
and keeps in the same state after the colhsion, the scat- 
tering is elastic and the scattering process associated to 
the ground state of the QD is of dominant contribution 
to the scattering cross-section. In this case, one-channel 
treatment can be a reasonably good approximation to 
calculate scattering cross-section even if the incident elec- 
tron is of enough energy and thus several inelastic chan- 
nels are open during the collision. When only the elastic 
channel is considered {i.e. i = m = n ^ 0), Eq. Q is 
reduced to 



2m* 



(v2-Kfc')V(r) = ^(r)V'(r), 



(18) 



where ^p{r) — V'oo(r), V{r) — Vbo(r) and k = ko. The 
L-S equation for the scattered electron becomes 

V'(r) = ^k(r) + / dr'G(°)(k,r,r')T^(r')^(r'), (19) 



where </?k(r) — e'*^^ and the Green's function is given by 

Eq. nni). 

When the scattering potential in the above equation 
is central, i.e., V{r) = V{r), the L-S equation can be 
solved easily using the partial wave expansion technique 
as described in Appendix B. Moreover, the scattering am- 
plitude or the cross-section in this case can be obtained 
in terms of the phase-shifts of the partial waves. The 
DCS as a function of partial wave phase-shift A; is given 
as: 



2 

nk 



Kie'^'^' sinA, cos(/6') 



and the ICS is given by 



Too (-Bo) ^ sin^ 

1=0 



where = 1 for Z = and ki — 2 ior I ^ 0. 



B. Scattering by a one-electron QD 

Electron scattering and transport through a QD of a 
few electrons are currently of great experimental and the- 
oretical interest. Here we present the case of a QD with 
only one confined electron. We focus on the exchange 
effect on the electron scattering and the spin-flip scat- 
tering mechanism. The total Hamiltonian Eq. ([2]) in the 
one-electron QD case is given by. 



ff(ri,r2) = 



2m 



VQD(r2) +i/QD(ri) +Mnt(ri,r2), 

(20) 



Mnt(ri,r2) 



£5 11-2 -ri I 



(21) 



where ri labels the localized electron in the QD and r2 
the incident electron. As we have mentioned, to solve the 
scattering problem, we need firstly to know the electron 
states in the QD which are determined by the following 
equation, 



i?QD(r)C"(r) 



2m 



:V2-fyQD(r) C"(r)=e„C"(r). 

(22) 

The solution of this one-electron QD is straightforward 
as soon as the confinement potential Vqd is defined. 

According to Eq. ([9]), there is an infinite number of 
quantum states involved in the scattering. In performing 
a numerical calculation, however, we have to truncate 
this to a finite number of states. As a matter of fact, 
when the QD is initially in its ground state and the in- 
cident electron has a small kinetic energy, it is a good 
approximation to consider only a few scattering chan- 
nels associated to the low-energy levels of the QD. In the 
present calculation, we consider the channels associated 
to the ground state eg and two excited states ei and £2 of 
the QD. When the incident electron passes through the 
QD initially in the ground state £0, the scattering can be 
either elastic keeping the QD in the same state or inelas- 
tic leaving behind the QD in an excited state. According 
to Eq. ([5]), the spatial part of the total wave-function of 
the two electrons (one incident and the other confined) 
can be written as, within a three-level model, 

2 

vl',;(ri,r2) = [C"(ri)V'™(r2) ± C (r2)V'™(ri)] , (23) 

n=0 

where the signs {+) and (— ) represent the singlet and 
triplet total-spin states of the two-electron system, re- 
spectively. The scattering equation [Eq. ([9])] becomes. 



2m^ 



(V 



V'mi(l"2) 



2 

E 

n=0 



[C„(r2)±yr„(r2)]V^n.(r2) 



with 



and 



c*„ = v^QD<5™„ + (nvi„tin 



(24) 



(25) 



21.2 



n^k: 



2m' 



/ rr) '■' 

where i and m {i,m = 0,1, and 2) indicate the initial 
and final state of the system, respectively. 

According to conservation of the total energy of the 
system, the relation between the kinetic energies of the 
incident (scattered) electron and the energies of the QD 
is given as 



£0 



2m* 



£1 



2m* 



£2 



n'kl 

2m* ■ 



(27) 
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TABLE I: Tangent of the phase shift of different partial waves I = 0, 1, 2, 3, and 4 for the first six iterations within the MCF for an 
incident electron of kinetic energy £"0=0.6 meV. The number between bracket represents the power of ten, e.g., (—4) = 10^^. 



iteration 





1 


2 


3 


4 


5 


6 


tan Ao 


-7.2103 


-6.0844 


-1.7589 


-1.6348 


-1.6442 


-1.6347 


-1.6347 


tan Ai 


-0.8833 


-0.6862 


-0.6874 


-0.6305 


-0.6305 


-0.6305 


-0.6305 


tan A2 


-0.0288 


0.1838 


0.2465 


0.2495 


0.2495 


0.2495 


0.2495 


tan A3 


5.43(-4) 


-2.62(-4) 


-3.66(-5) 


-3.64(-5) 


-3.64(-5) 


-3.64(-5) 


-3.64(-5) 


tan A4 


1.23(-4) 


-1.20(-4) 


-1.20(-4) 


-1.20(-4) 


-1.20(-4) 


-1.20(-4) 


-1.20(-4) 



The corresponding L-S equation is reduced to 

+ /dr'Go(k™,r,r')F„„(r')V'™(r'), (28) 

where Vmn = ^^mn i ^nm- The different channels are 
coupled through the potential matrix elements Vmn of 
the same total energy. In other words, the scattering for 
an incident electron of momentum fcg (associated to the 
QD ground state eo) couples to that for an incident elec- 
tron of momentum ki (associated to the first excited QD 
state £1) satisfying Eq. ([27ll . Because the two electrons 
can form both the singlet (+) and triplet (— ) states, the 
scattering cross-sections are different for these two dis- 
tinct cases: 

-™(^) = ^l/i!k(^)l' (29) 
for the singlet state, and 

^Uo) - '^\.fi\M' (30) 

for the triplet states, where the scattering amplitudes are 
given by 

x[y^*„(r')±C^„(r')]^™(r'). (31) 

The total differential cross-section or the spin- 
unpolarized (su) DCS is determined by a statistical ad- 
mixture of the singlet and triplet state scattering, 

<"(^) = jKn.(^)+3aL(0)), (32) 

where the factor 3 in the equation is due to statistical 
weight of triplet states. Another interesting quantity is 
the spin-flip (sf) DCS which describes the spin-flip scat- 
tering probability of an incident electron resulting from 
the exchange interactior>i^. The sf-DCS is given by. 



where 

/feL,fe.(^) = /l!k(^)-/l:k(f^) 

= -\\[^il I dr'e-'*'" "-V-(r')V™(r'). (34) 

IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

We model the conflnement potential of the QD by a 
2D finite parabolic potential 

^QD(r) = | ' (35) 

\ , r > ro, 

where wq is the conflnement frequency and tq is the ra- 
dius (size) of the dot. We will calculate in this section 
the scattering due to a one-electron quantum dot. For 
such a system, the solution of Eq. is straightfor- 
ward. We expand the eigenfunction ^" in the Fock- 
Darwin basis22i and diagonalize numerically the Hamil- 
tonian. The eigenstates can be labeled by a set of quan- 
tum numbers n—{j,l) with the radial quantum number 
j=0,l,... and the angular momentum quantum number 
1=0, ±1,... . The state (j = 0,1 = 0) is the ground state 
n = and the flrst two excited states {j = 0,1 = ±1) are 
degenerate corresponding to n = 1 and n = 2 [ei = £2). 

In order to solve the L-S equations, we use the par- 
tial wave expansion in two dimensional system combined 
with the MCF. All the involved functions are expanded 
in the angular momentum basis so that we obtain a radial 
L-S equation for each angular momentum. Numerically, 
we are able to choose the components of the angular mo- 
mentum which contribute to the cross-sections up to a 
desirable precision. In Appendix B we show how the par- 
tial wave expansion can be applied to the multichannel 
L-S equations in a two-dimensional system. 

A. Convergency of the MCF 

The MCF was applied to the electron-atom 
scatteringi^ and electron-molecule scattering!^. In 
all those cases, it has shown a rapid convergency. Here, 
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FIG. 1: (Color online) The elastic DCS's obtained within the one-channel model for electron scattering by the one-electron QD of Hluq = 5 
meV and ro = 35 nm. The incident electron energies are indicated in the figures, (a) The spin-unpolarized DCS with (the solid curves) 
and without (the dashed curves) the exchange potential; (b) The DCS due to the singlet state (the solid curves) and the triplet state (the 
dashed curves); and (c) The spin-flip DCS. 



we apply the MCF to the electron-QD scattermg in a 
two-dimensional semiconductor nanostructure. First of 
all, we check the convergency of this method for electron 
scattering through a QD. We consider a one-electron QD 
with hujo = 5 meV, tq = 35 nm and an incident electron 
of kinetic energy Eq = h^kQ/2m* = 0.6 meV. The results 
are obtained within the one-channel approximation. For 
simplicity, only the static scattering potential is consid- 
ered and the exchange potential is neglected in Eq. ((3T|l . 
Table I gives the calculated partial wave phase-shifts, for 
angular momenta up to Z = 4, of the first six iterations. 
Because the localization length of the con fined elec tron 
wave-function in the QD is about uq = •\/?i/m*wo, the 
screening parameter is taken as X — a^^ throughout this 
paper, (ag = 14.75 nm for a GaAs QD of Tiluq — 5 meV). 
We have performed calculations with different values of 
A. The calculated results have showed that, although a 
smaller A has enhanced the static potential scattering, 
the exchange effects and the spin-flip scattering are not 
affected significantly. From Table I we can see that the 
phase-shifts converge at the fifth iteration. It also shows 
that the first Born approximation, which corresponds 
to our zero-th iteration calculation, is indeed a very 
poor approximation in dealing with the electron-QD 
scattering. In order to obtain a correct scattering 
cross-section through a QD, it is necessary to use a 
robust method such as the MCF. Although the results in 
Table I are for a particular case, we have verified that all 
our calculations (with or without exchange interaction) 
in the present paper are convergent within 6 iterations. 



B. One-channel scattering 

Within the one-channel approximation, we calculated 
the elastic DCS for electron scattering by the one-electron 
QD of tkuQ — 5 meV and tq =35 nm at incident kinetic 



energies Eg =0.6, 1.7, and 4.2 meV. The obtained DCS's 
are presented in Fig. 1 as a function of the scattering 
angle. Fig. 1(a) shows the total or the su-DCS (Jqo{0). 
To illustrate the effect of exchange interaction, the su- 
DCS's due to the static potential only (i.e., neglecting 
the exchange potential in Eqs. and ([?!]) ) are given 
by the dashed curves in the figure. We see that the ex- 
change interaction is of significant contribution to the 
low-energy and/or small angle scattering. The exchange 
effect on the scattering is originated from the two dif- 
ferent coupling states between the incident and the QD 
electrons (i.e., the singlet and the triplet states) during 
the collision as indicated in Eq. ((24l) . The correspond- 
ing DCS's due to the singlet [crQQ((?)] and triplet states 
[(Tqq{9)] defined by Eqs. (^5]) and ([50]) . respectively, are 
shown in Fig. 1(b). In Fig. 1(c), we plot the spin-fiip 
(sf) DCS crf){){9) given by Eq. Comparing Fig. 1(a) 

with Fig. 1(b), one can see that the exchange interaction 
affects more strongly the su-DCS when the difference be- 
tween <Tqq{9) and (Toq{0) is large. In Fig. 1(c) we observe 
that the spin-ffip scattering due to exchange potential 
occurs mostly around 9 ~ 90° at lower incident energy 
{Ea =0.6 meV). For higher energies {E^ =1.7 and 4.2 
meV), the spin-fiip scattering is more relevant at small 
scattering angles. 

Fig. 2(a) shows the ICS as a function of incident elec- 
tron energy Eq for the spin-unpolarized scattering and for 
that considering the static potential only. We see that the 
exchange interaction affects significantly the ICS at low 
Eq. At higher energies, however, the ICS is dominated 
by the static potential. In Fig. 2(a) we also present the 
spin- flip ICS (the dashed curve). A maximum spin-fiip 
probability is found at Eq =1.1 meV which is about 37% 
of the total scattering. In Fig. 2(b), we plot the ICS due 
to the singlet and the triplet states. It shows a strong 
dependence of the ICS on the spin states of two electrons 
in the system. 
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FIG. 2: (Color online) The elastic ICS as a function of Eq for the one-electron QD. (a) The su-ICS (the solid curve), the sf-ICS (the 
dashed curve), and the ICS due to the static potential only (the dotted-dash curve); (b) The ICS due to the singlet (the solid curve) and 
triplet states (the dashed curve) during the scattering. 




FIG. 3: (Color online) The phase-shift as a function of Eq for the partial waves / = 0, 1, 2, 3, and 4 due to the static potential scattering 
of the one-electron QD. 



The scattering peaks in ICS due to the static potential 
(the dotted-dash curve in Fig. 2(a)) at Eq — 1.22 and 
6.0 meV are due to the occurrence of the so-called shape 
resonances, resulting from a virtual confined state at the 
corresponding energy. In order to clarify the origin of 
these features, we plot in Fig. 3 the corresponding partial 
wave phase-shifts A/ (for ^ = 0, 1, 2, 3 and 4) due to the 
static potential. At Eq — > 0, Aq and Ai are larger than 
7r/2 indicating the presence of the bound states of angular 
momenta I =0 and 1 in the QD. A rapid increase of A2 
(A3) at around Eq = 1.2 mcV (6.0 meV) corresponds to 
a virtual bound state of ^ = 2 (^ = 3) leading to the shape 
resonance scattering peak in the ICS. Similarly, peaks in 
the ICS at £;o = 1.72 meV (0.57 meV) for the singlet 
(triplet) state scattering in Fig. 2(b) can be related to 
the virtual bound states in the system. The broad peak 
in the ICS of triplet state is possibly result from a virtual 
state of two interacting electrons. 



C. Multi-channel scattering 



The energy difference between the first excited state 
and the ground state is 4.90 meV in the QD of hLLjQ=5 
meV and ro=35 nm. When the kinetic energy Eq of an 
incident electron is higher than this energy difference, the 
inelastic scattering channel is open which leaves the QD 
in the excited state £1 after scattering. In such a case, 
the multichannel scattering process has to be considered. 
When the three lowest states of a one-electron QD are 
included in the calculation, there are 9 possible scatter- 
ing channels. For the present QD, as the first excited 
state is two- fold degenerate, i.e. Si —62, we find the fol- 
lowing scattering cross-sections: the elastic roo(i?o) and 
inelastic scattering rio(£'o) for the QD initially in its 
ground state; the elastic scattering rii(ii^i) and super- 
elastic scattering roi(-Ei) for the QD initially in the first 
excited state. There are also two excitation cross-sections 
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FIG. 4; (Color online) The multichannel integral cross-sections (the solid curves) and T'^^ (the dashed curves) as a function of Ei 
[i = 0, 1) for the one-electron QD considering the three lowest energy states. The thin curves in (a) are the corresponding results within 
the one-channel model. 



^2i{Ei) — Ti2{Ei) for the QD in the excited states of 
different angular momenta, although in these cases the 
energy of the scattered electron is the same as that of 
the elastic scattering rii(_Ei), i.e. e\ — £2 and fci — k-i 
in Eq ([27|) . In Fig. (4) we show the different integral 
cross-sections due to the singlet and the triplet states. 
Figs. 4(a) and 4(b) give the elastic ICS (Foo) and in- 
elastic ICS (Fio), respectively, for the QD initially in its 
ground state. The inelastic cross-section is two orders 
of magnitude smaller than the elastic one. Moreover, the 
inelastic scattering due to the triplet state is very weak at 
higher energies. Coupling between different QD levels (or 
different scattering channels) leads to resonant scattering 
on both the elastic and inelastic cross-sections. The thin 
curves in Fig. 4(a) are the previous results in Fig. 2(b) of 
the ICS within the one-channel model. We see that the 
one-channel approximation yields quite good results for 
the elastic scattering. 



Figs. 4(c) and 4(d) show the elastic ICS for the QD in 
the first excited state. At small incident energy £'1, the 
scattering due to the triplet state is much stronger than 
that due to the singlet state in this case. If the QD re- 
mains in the same excited state after the scattering, the 
ICS Fii (=r22) as shown in Fig. 4(c) is large at higher en- 
ergy and decreases slowly with increasing E\. However, 
if the angular momentum changes after the scattering, 
the ICS r2i (=ri2) vanishes rapidly (both due to the 
singlet and triplet states) at high incident energies. 



V. CONCLUSION REMARKS 

We presented a general method to calculate the elec- 
tron scattering through an N-electron QD embedded in 
a two-dimensional semiconductor system. The multi- 
channel L-S equations are solved numerically using the 
iterative method of continued fractions considering the 
electron-electron interactions. We applied this method 
to the case where only one electron is confined in the QD. 
The results indicate a rapid convergency of this method 
for two-dimensional scattering in a semiconductor nanos- 
tructurc. It also shows that the first Born approxima- 
tion is so poor that cannot yield correct scattering cross- 
section. 

We found that the exchange effects are relevant when 
the kinetic energy of incident electron is small, as showed 
by the obtained DCS and ICS. The shape resonances were 
found in the elastic ICS including or not the exchange po- 
tential. The spin-flip cross-section due to exchange inter- 
action shows a maximum both in the DCS as a function 
of the scattering angle and in the ICS as a function of the 
incident electron energy. The maximum spin-flip scatter- 
ing reaches as high as more than 30% in comparison to 
the total scattering. In multichannel scattering including 
the excited states of the QD, we obtained the inelastic 
scattering cross-sections. They are about two orders of 
magnitude smaller than the elastic ones. 

In this paper, we emphasize the theoretical approach 
and numerical method to calculate the electron scattering 
by a charged quantum dot. The scattering cross-sections 
were obtained for a spin unpolarized system. It can be 
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extended to a spin polarized system which is of great 
interest for electron transport in semiconductor nanos- 
tructures. The application to a spin polarized system is 
straightforward as soon as the initial spin states of the 
system are defined. In the numerical calculation, we pre- 
sented the cross-sections due to a one-electron QD scat- 
tering. For a QD of two or more electrons, we need firstly 
to know the eigenstates of the QD with electron-electron 
interactions. Then the scattering cross-sections can be 
calculated according to the total wavefunction defined 
by Eq. ^ as what has been done for electron-atom and 
electron-molecule scattering where several electrons are 
presented 
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APPENDIX B: PARTIAL WAVE EXPANSION 

In two dimensions the angular momentum basis is 
given by2^, 



(Bl) 



where I = 0,1,2,..., ki — 2 for / 7^ and ki — 1 for 
1 = 0. In applying the partial wave expansion in the 
multi-channel scattering problem Eq. (jl2p . we expand 
all functions, i.e., the incident free electron wavefunc- 
tion the Green's function G''-°-' (km, r, r'), and the 
scattered electron wavefunction '0mi(r), in the angular 
momentum basis as follows. 



ft 



00 I— 

W= E J^^'Mkr)6u'ei{c^r)ei'{<l>k), (B2) 



and 



APPENDIX A: METHOD OF CONTINUED 
FRACTIONS 



^™(r) = ^ V'^';(fc,r)e,(0,)eH0fc), (B3) 



Ll'=0 



The method of continued fractions (MCF)ii is an it- 
erative method to solve the L-S equation. To apply this 
method for a multi-channel scattering we have firstly to 
rewrite Eq. I|12p in a matrix form: 



f- 



(Al) 



In the first step to start the MCE, we use the scattering 
potential V = V^^"^ and the free electron wave-function 
in Eq. (jAl[) . Afterwards, we define the nth- 



\f 



(o)\ 



order weakened potential as 



where 



(A2) 



(A3) 



The ni/i-order correction of the T matrix can be obtained 
through 

(^(") |^(»)) _ ((^(") |^(")) . (A4) 

Hence, we can stop the iteration when the potential y(^) 
becomes weaker enough. In the numerical calculation, we 
start with T^^+i) = and evaluate t'^n) ^j,(n -1) ^ ^^^^ 
T^^\ Therefore the T matrix is given by 

(A5) 



where 0^ and (j)k are the variables due to expansion on the 
position r and momentum k, respectively. The expansion 
on the Green's function yields the following expression. 



G(")(k,„,r,r') = 



where k = ki, r< = min(r, r'), r> = max(r, r'), 
Ji{kmr) {Yiikmr)) is the Bessel (Neumann) function 
and Hi^\kmr) = Ji{km^) + iYi{kmr) is the Han- 
kel function^^. Using the partial wave expansion the 
Lippmann-Schwinger equation can be reduced to a set of 
radial equations. The radial Lippmann-Schwinger equa- 
tion corresponding to Eq. (|12p is given by. 



(B4) 



&,r) = J^t'Mkr)6u'6^ 



(B5) 



00 00 „oo 

EE / r'dr'gl{k^,r,r')V:^':{r')^2\r'), 



where 



g^(fcm,r,r') = _ J-1 J,(fc„r<)i?/i)(fcmr>) (B6) 



and 



Vhn{r')= I dcj,r'Ql{<Pr')V^n{v')Qv'{^r'). (B7) 

Jo 

We see that, when the partial wave method is used, there 
is a change in the continuum variable to a partial wave 
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I. Consequently, the wavefunction ipmii'<^) becomes a ma- 
trix function with elements tpmii^jr). 

The partial wave expansion for the exchange potential 
is a little subtle due to its non-locality. Here, we show 
some details how the partial wave expansion is applied in 
this case. We take as an example the exchange potential 
which couples the channels n and m for a single electron 
spin-orbital a, 

WV-n.W = -4CW / driC*(ri)T-^^™(ri). 

|r ri| 

(B8) 

The partial wave expansion of the spin-orbital function 
is given by 



(B9) 



(=0 



The product of two different functions can also be ex- 
panded in the angular momentum basis as follows, 

V'^wcrw ^Y.^^no.irm'f>r)Qi'{M, (biq) 



where 



rf''-' (r) 



E 

A, A' 



A+A' 



(Bll) 



Using the above relation, we obtain Eq. (jBSp in the par- 
tial wave expansion form, 

F-WV',^^ = (B12) 

1,1' 

X / ndnU^^ in) / _========, 

Jo Jo \Jt'^ + — zrri cos^W) 

where 9 — <j)r — 4>ri ■ To solve the angular integral we use 
the generating function of the Legendre Polynomials^^, 



obtain finally the exchange potential 

^0 1,1' 
°° foo , j 

xJ2 ndnn'^.^^^{n)ci,,-^. (Bi5) 

In the numerical calculations, wc firstly evaluate the co- 
efficients ci j given by Eq. (|B14[) . Then the integration 
on ri in Eq. (jB15[) is performed for each iteration in the 
MCF. Finally we multiply the result by -f^C2(r)- 

Within the one-channel approximation (i = m = n = 
0) , the calculations can be further simplified by using the 
concept of phase shift. Considering a central potential 
V{r) {l^V = I"), Eq. dm becomes 



^'(fc,r) = ^l^iiJi{kr)+ I r'dr'gly{k,r,r')V{r')tl;\k,r') 

(B16) 

where il)^{k^r) = ipQQ^kjr). To define the phase-shift we 
write the asymptotic form of the above equation as 



V;'(fc,r) ^ Ai^l.cos{kr-^-^~^-Ai), (B17) 



where A/ is the phase-shift. Comparing the coefficients 
of eJ'^^ and e"*'^'' of Eq. (|B17|) with the asymptotic form 
of Eq. (jB16p one can obtain the following relations 



(B18) 



and 

e'^' sin A, 



— TT 
27 



r'dr'Ji{kr')V{r')i;\r'). (B19) 



On the other hand, from the definition of the scatter- 
ing amplitude in Eq. (I14p . we can express the scattering 
amplitude fko,ko in terms of the phase-shift^^ A/, 



1 



^j,2 _|_ J.2 „ 2rri cos{9) 



j- 

3=0 ^> 



COS ( 



(B13) 



where r< = min(r, ri), r> = max(r, ri) and Pj {cos 0) 
are the Legendre Polynomials. Thus the angular integral 
that we need to solve is 



ho,koiO) = 2 V sin A,e,(^) 



1=0 



(B20) 



The corresponding DCS is croo(^) = l/fco,fco(^)l ^^d 
the ICS is given by 



deQi{9)Pj{cost 



(B14) 



Substituting the Eqs. (|BT3)) and ([BT4)) into Eq. (|BT2)) we 



4 2 a 

Poo = T 2^ sm A;. 



(B21) 



1=0 
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